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Abstract 

Mesophyll conductance (g m ) can markedly limit photosynthetic C0 2 assimilation and is required to estimate the param- 
eters of the Farquhar-von Caemmerer-Berry (FvCB) model properly. The variable J (electron transport rate) is the 
most frequently used method for estimating g m , and the correct determination of J is one of its requirements. Recent 
evidence has shown that calibrating J can lead to some errors in estimating g m , but to what extent the parameteriza- 
tion of the FvCB model is affected by calibrations is not well known. In addition to determining the FvCB parameters, 
variants of the J calibration method were tested to address whether varying C0 2 or light levels, possible alternative 
electron sinks, or contrasting leaf structural properties might play a role in determining differences in ap, the product 
of the leaf absorptance (a) and the photosystem II optical cross-section (P). It was shown that differences in ap were 
mainly attributed to the use of A/C, or A/PPFD curves to calibrate J. The different ap values greatly influenced g m , lead- 
ing to a high number of unrealistic values in addition to affecting the estimates of the FvCB model parameters. A new 
approach was devised to retrieve leaf respiration in the light from combined A/C, and A/C c curves and a framework to 
understand the high variation in observed g m values. Overall, a background is provided to decrease the noise in g m , 
facilitating data reporting and allowing better retrieval of the information presented in A/C, and A/C c curves. 

Key words: A/C, curve fitting, chlorophyll fluorescence, Coffea arabica, Limonium gibertii, Nicotiana tabacum, variable J 
method. 



Introduction 

Photosynthesis is a major process that affects plant growth 
and crop productivity. In addition to stomatal and biochemi- 
cal factors, the photosynthetic capacity of leaves is also deter- 
mined by the mesophyll conductance (gm), which regulates 
the C0 2 flux from the intercellular airspaces to the sites of 
carboxylation in the chloroplastic stroma (Flexas et al, 2012). 
Early gas exchange studies assumed infinite and constant g m , 
implying that the C0 2 concentrations in the substomatal cavi- 
ties (Q) and chloroplasts (C c ) would be the same (Farquhar 
et cil., 1980). However, the role of finite g m limiting photosyn- 
thesis in response to several biotic and abiotic stresses (Flexas 



et al., 2008), as well as its importance in constraining maxi- 
mum photosynthetic rates, particularly in evergreen sclero- 
phylls (Warren et al. , 2004; Niinemets et al. , 201 1), is now well 
established. Actually, there is a general consensus that g m must 
be incorporated into the Farquhar-von Caemmerer-Berry 
(FvCB) model of leaf photosynthesis (Farquhar et al., 1980) 
because this model underestimates the maximum Rubisco 
carboxylation rate ( K cmax ) when g m is considered to be infinite 
(Ethier and Livingston, 2004; Niinemets et al., 2009a). 

Several methods to estimate g m have been reported 
(Warren, 2006; Pons et al., 2009). Among them, those based 
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on gas exchange coupled with the chlorophyll fluorescence 
technique have been extensively used. In particular, the vari- 
able / method (where / stands for electron transport rate, 
usually measured using chlorophyll fluorescence analysis) 
is the most widespread method because it allows the point- 
specific estimation of g m as well as tracking how g m suppos- 
edly changes in response to varying light or C0 2 (Pons et al. , 
2009). However, the variable / method is highly sensitive to 
errors in some parameters, especially the C0 2 compensation 
point in the absence of photorespiration (T*) and / (Harley 
et al, 1992). Importantly, / needs to be calibrated due to the 
uncertainties in both the leaf absorptance (a) and photosys- 
tem II (PSII) optical cross-section ((3). 

The /calibration is based on the analysis of response curves 
of photosynthetic rates to light L4/photosynthetic photon 
flux density (PPFD)] or C0 2 (A/Q) under non-photorespi- 
ratory conditions, typically under low (1-2%) oxygen condi- 
tions. Under these circumstances, the relationship between / 
calculated through gas exchange (/a) and chlorophyll fluo- 
rescence (/ F ) is expected to be linear because electron trans- 
port flow is primarily associated with Rubisco carboxylation. 
Consequently, J A can be calculated as J A -4(A+R L ) (Warren, 
2006). Other assumptions are that a(3 or light respiration 
(R L ) values do not change in response to varying Q or PPFD. 
To date, three different approaches have been used to cali- 
brate /: the relationship between the quantum yield of PSII 
(0 PSII ) and C0 2 (0 C02 ) (Valentini et al, 1995), the relation- 
ship between / A and / F (Pons et al , 2009), and the plot of 
A versus PPFD O psii /4, which is also used to estimate R L 
(Yin et al, 2009). Another point deserving attention is the 
use of AfPPFD or Al Q curves under low 0 2 to perform the 
calibration; ideally, both curves should give the same results. 
However, as observed by Yin et al. (2009), fundamental dif- 
ferences exist when choosing ^4/PPFD or A/Q curves to 
calibrate /. First, when using the full A/Q curve, photores- 
piration can still occur at low C0 2 levels. Alternatively, the 
excess of energy can increase the rate of alternative electron 
flow. Both cases can compromise the linearity of the relation- 
ship between J A and / F . The same problem occurs in /4/PPFD 
curves, where higher PPFD intensities can induce alternative 
electron flow (depending on g s , low Q can occur as well). To 
overcome this problem, Yin et al (2009) recommended using 
the electron transport-limited regions of both curves, namely 
the combination of low light levels from ^4/PPFD and high Q 
levels from Al Q, which would ultimately minimize the risk of 
photorespiration and alternative electron flow. 

More recently, Gilbert et al. (2012) examined the use of 
^4/PPFD or Aid curves to calibrate / and demonstrated 
dramatic changes in g m values estimated using either the Al 
PPFD or A/Ci curves. However, to what extent these cali- 
brations affect estimates of photosynthetic parameters such 
as J^cmax or / max , and whether this effect is dependent on 
the species studied, remain unresolved. Flexas et al. (2007) 
found no differences using both calibrations for tobacco, 
and extrapolated this result to other species where only the 
^4/PPFD calibration was used. Hassiotou et al. (2009), using 
the calibration based on A/Q curves performed under two 
light levels, concluded that calibration is light dependent. 



Other studies did not use the ^4/PPFD or A/Q calibration, 
but estimated a using integrating spheres and assumed a (3 
value of 0.5 (e.g. Galmes et al, 2007; Tosens et al, 2012), 
even though (3 has been reported to vary (Laisk and Loreto, 
1996). Collectively, these results reveal no consensus on how 
to calibrate /, even though g m estimated by the variable / 
method is mainly dependent on /. In addition, because the 
fluorescence signal primarily emanates from the upper meso- 
phyll layers, but gas exchange parameters are volume based 
(Warren et al, 2006), measured / and estimated g m are most 
probably less representative of the whole leaf in species with 
lower specific leaf area (SLA; leaf area per unit dry mass). 

The main questions asked in this study were as follows, (i) 
How do the / calibrations affect g m l (ii) How might these cali- 
brations be translated into the FvCB parameters? (iii) Would 
the calibrations be dependent on species? To answer these ques- 
tions, measurements were carried out using three species with 
contrasting SLA and photosynthetic capacities. The results 
highlight how g m , F cmax , and / max may be affected by using Al Q 
or ^4/PPFD curves under low 0 2 to calibrate /. The results are 
discussed in the context of current models of g m estimations. 



Materials and methods 

Plant material and growth conditions 

Limonium gibertii (Senn.) Senn. and Nicotimw tabacum L. seeds 
were germinated, and plants were grown outdoors under typical 
Mediterranean climate conditions (Balearic Islands, Spain, 39°38'N, 
2°38'E, 85 m a.s.l.) in 4 litre pots with a commercial substrate (horti- 
cultural peat) and perlite at a proportion of 4:1. Seedlings of Coffea 
arabica L. obtained from seeds were grown outdoors under subtrop- 
ical conditions in Vicosa (20°45'S, 42°15'W, 650 m a.s.l), southeast- 
ern Brazil, using 12 litre pots containing a mixture of soil, sand, 
and composted manure (4:1:1, v/v/v). Plants were irrigated and ferti- 
lized as required. Measurements were performed during the summer 
(growing season) on ~l-year-old plants in the case of C. arabica and 
L. gibertii, and 1 -month-old plants in the case of N. tabacum, on 4-6 
plants per species. 



Gas exchange and fluorescence measurements 

Leaf gas exchange and chlorophyll a fluorescence were measured 
simultaneously with an open-flow infrared gas-exchange analyser 
system equipped with a leaf chamber fluorometer (LI-6400XT, 
Li-Cor, Lincoln, NE, USA). Environmental conditions in the leaf 
chamber consisted of a leaf-to-air vapour pressure deficit of 1.2-2.0 
kPa and a leaf temperature of 25 °C. 

In light-adapted leaves, the actual O psii was determined by meas- 
uring steady-state fluorescence (F s ) and maximum fluorescence dur- 
ing a light-saturating pulse of ~8000 umol nr 2 s~' (F m "), following 
the procedures of Genty et al. (1989): 

0>P Sn =(/' m '-F s )/^ m ' (1) 

The electron transport rate (J F ) was then calculated as: 

7 F =ap PPFD O psii (2) 

where PPFD is the photosynthetically active photon flux density, 
a is the leaf absorptance, and |3 is the PSII optical cross-section. The 
product a|3 was determined from the relationship between O psii and 
flfco, or ^ an d PPFDOpsn/4, obtained by varying either the light 



Improving the estimation of mesophyll conductance | 3287 



intensity or the C0 2 concentration under non-photorespiratory con- 
ditions in an atmosphere containing <1% 0 2 (Valentini et al, 1995; 
Yiner al., 2009). 

Four to six A/Q and ^/PPFD curves under <1% 0 2 (,4/PPFD and 
A/Q) or 21% 0 2 (only A/Q) were obtained from different plants for 
each species. In light-adapted leaves, A/Q curves were initiated at an 
ambient C0 2 concentration (Q) of 400 umol moF 1 under a saturating 
PPFD of 1500 umol nr 2 s -1 . Once steady state was reached, C a was 
decreased stepwise down to 50 umol moF 1 air. Upon completion of 
the measurements at low C a , it was returned to 400 umol moF 1 air to 
restore the original A. Next, C a was increased stepwise to 2000 umol 
moF 1 air. For the ^4/PPFD curves, C a was held at 400 umol moF 1 , and 
the curve was initiated at a PPFD of 1500 umol nr 2 s _1 ; then, PPFD 
levels were decreased to 0 umol nr 2 s . Both the A/Q and ^/PPFD 
curves consisted of 11-13 different C a values or PPFD intensities. 

C c was calculated after Harley et al. (1992) as: 



In addition, the method based on Yin's approach was also used 
(Yin et al. , 2009), which presents a straightforward way to derive a|3 
and R L , as follows: 

J A =4 (A+R L ) (8) 
/ F =ap PPFD 0 PSII (9) 

J A is rewritten as 

A=J;J4~R L (10) 
Assuming J A -J F under non-photorespiratory conditions gives 

A=a$ PPFD Opsn/4-^L (11) 



C c ={r*[/ F +8(,4+ J R L )]}/[/ F -^t+ J R L )] (3) 

where Y* was determined from the in vitro Rubisco specificity fac- 
tor (S do ) (see below) as: 

r*=o/s c/0 (4) 

A was taken from gas-exchange measurements, and the J F val- 
ues were obtained from chlorophyll a fluorescence yield. The rate 
of mitochondrial respiration at darkness (-R dark ) was measured early 
in the morning in dark-adapted leaves, and it was divided by two 
(-Rdark/2) to serve as a proxy for R L . 

After estimating Q, g m was calculated as follows (Harley et al. , 
1992): 

g m =AI(Q-Q) (5) 

From the A/Q and AIC C curves, the maximum carboxylation 
capacity (K cmax ) and maximum capacity for electron transport rate 
(/ max ) were calculated on a Q and C c basis using the kinetic param- 
eters of Rubisco described below and, for comparative purposes, 
those described in Bernacchi et al. (2002). The FvCB model was fit- 
ted to the data by applying iterative curve fitting (minimum least 
square difference) using the Microsoft Excel Solver tool (Microsoft 
Corporation, Redmond, WA, USA). Additionally, g m , ^ and 
/ max were estimated using the Ethier and Livingston (2004) method, 
which is based on fitting Al Q curves with a non-rectangular hyper- 
bola version of the FvCB model, relying on the hypothesis that g m 
reduces the curvature of the Rubisco-limited portion of an A/Q 
response curve. For the method based on fitting A/Q curves, species- 
specific S do values were used as in the Harley method. Corrections 
for the leakage of C0 2 and water vapour into and out of the leaf 
chamber of the Li-6400^40 have been applied to all gas-exchange 
data, as described by Rodeghiero et al. (2007). The percentage cor- 
rections applied to C0 2 and water vapour flux rates are shown in 
Supplementary Table SI available at JXB online for the different 
species. 

Calibration relationships 

The relationship between J A and J F was calibrated using the linear 
plot of <t> Psn and <1> C02 , based on Valentini et al. (1995): 

3l» SII = k% 02 +b (6) 

/ F =4(0) PSIr -b) PPFD/Ar (7) 

where 4/k=afi. k and b were obtained through the linear fit of 
0 PSII versus <l>co 2 ■ 



As Equation 1 1 has the form of y-ax+b, through the linear fit of 
A versus PPFDO PSII /4, R L can be retrieved as the j'-intercept, and 
a|3 can be retrieved as the slope of the regression. Notably, this equa- 
tion presented by Yin et al. (2009) has an extension of the FvCB 
model to account for alternative electron fluxes in the form of pseu- 
docyclic (/ pseu d 0 ) and cyclic (/ cyc ) electron flow (see further details in 
Yin et al., 2004): 

A=a$ PPFD (D PSn { l-[/p Seu do/(l-/c y c)]}/^L (12) 

According to the updated model, to accomplish Equation 8, 
not only are non-photorespiratory conditions required but also 
the down-regulation of alternative electron fluxes so that low val- 
ues of /p S eudo and / cyc allow Equation 12 to be as close as possible to 
Equation 11, thus reliably estimating a|3. 

Four types of calibration were devised based on the approaches 
of Valentini and Yin described above. The first and second calibra- 
tion methods were based on the Op^/ <&co 2 relationship, and the 
third and fourth methods were based on Yin et al. (2009) as follows: 



(i) Opsn/ <E>co7 (A/Cj): this calibration is performed using the entire 
A/Q curve under low 0 2 , and the PPFD level is that used in the 
normal A/Q curve (i.e. 1500 umol nr 2 s" 1 ). The R L value must 
be assumed. In this variant, the initial part of the A/Q curve 
is the most susceptible region for the occurrence of alternative 
electron sinks due to a high reductant (ATP and NADPH) sup- 
ply associated with high PPFD and a limitation on reductant 
use (low C0 2 and 0 2 ). Conversely, this variant most resembles 
the conditions used in the normal A/Q curve. 

(ii) O psii / © C(> (PPFD >400 umol nr 2 s -1 ): this is the variant more 
often reported in the literature and consists of using AIPPFD 
curves under low 0 2 and ambient C0 2 . A disadvantage of this 
variant is that a loss of linearity occurs at low PPFD, as © co 
is affected disproportionately by errors in respiration estima- 
tions or by the increase in mitochondrial respiration at low light 
(the Kok-effect; Brooks and Farquhar, 1985). Thus, it is recom- 
mended to exclude €> C07 values >0.05 to keep the relationship as 
linear as possible (Seaton and Walker, 1990; Edwards and Baker, 
1993). To meet this criterion, PPFD levels below 400 umol nr 2 
s" 1 had to be excluded, making this method prone to a high 
reductant supply (high PPFD) and intermediate reductant use 
(moderate C0 2 and low 0 2 ). 

(iii) Yin (PPFD <400 umol nr 2 s -1 ): this variant was originally 
described by Yin et al. (2009), and it can also be used to esti- 
mate R L . As this method is not based on quantum efficiency 
plots, there is not the disadvantage of having to exclude the 
points at low PPFD or assume a given R L value. Actually, it 
is recommended to use only the points at the linear phase of 
the A/VPFD (corresponding to PPFD levels <400 umol nT 2 
s~' for the species grown outdoors), which is in the range of a 
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lower reductant supply and where only the basal components 
of alternative electron sinks exist, 
(iv) Yin (PPFD <400 umol nr 2 s^/Q >500 umol moL 1 air): this 
method has the same advantages and assumptions of method 
(iii). A particularity of this approach resides in the incorporation 
of data from the A/Q curve with Q above 500 umol mor 1 air. 

Retrieving the respiration value from combined A/C, and 
A/C c curves 

The FvCB model predicts that A/Q curves and their A/C c counter- 
parts share the same C0 2 compensation point (r), as this parameter 
is independent of g m : 

f=[f*+k c ( 1 +o/k 0 )rj v^wx-Rd v cmax ) 

Evans and von Caemmerer (1996) proposed that an averaged g m 
can be estimated as the slope of the plot of A versus Q-C c . Because 
the plot goes through the origin, both A and Q-C c must be zero, 
which is what occurs when F is shared by both A/Q and AIC C curves, 
as predicted by the FvCB model. 

It was realized that for some data sets, when R L was set in advance 
as required to obtain C c , the plot of A versus Q-C c in the region 
which was Rubisco limited showed good linearity, although with 
an intercept differing from zero. Given the predictions of the FvCB 
model, this difference should be attributed to a biased R L used to 
calculate C c , which in turn leads to different F values, as obtained 
from the A/Q and A/Q curves. Given this fact and that the F* val- 
ues in this study were reliably estimated from Rubisco kinetics in 
purified preparations, it is proposed here that R L can be obtained 
as the value that makes the intercept of the plot of A versus Cj-C c 
equal zero. The approach was tested using ideal data sets where the 
respiration values used as inputs were successfully retrieved from 
the combined A/Q and A/Q curves. Consequently, for each curve, 
R L was estimated as described below. First, R L was set to be equal 
to -Rdark/2) as set in previous studies (Niinemets et al, 2005, 2006, 
2009ft). To maintain the plot A versus Q-Q as linear as possible, 
only points in the Q range strictly limited by Rubisco (Q <200 umol 
mol -1 air) were considered. Secondly, the new R L was obtained as 
the value that forces the intercept of the plot A versus Q-Q to be 
zero using the Goal Seek function available in Microsoft Excel (see 
Fig. 1 for an example). The new respiration value obtained in this 
way is hereafter referred to as R ACj /ACc ■ An Excel spreadsheet (with 
modelled and real data) is available (see Supplementary Spreadsheet 
SI at JXB online) that shows how this method was performed. 



Rubisco kinetic parameters 

The Rubisco kinetic parameters used in this study were measured 
in vitro, except the S c/o values for N. tabacum and L. gibertii, which 
were taken from Galmes et al. (2006) and Galmes et al. (2005), 
respectively. All measurements for the determination of Rubisco 
kinetic parameters were conducted at 25 °C The Rubisco speci- 
ficity factor was measured in highly purified extracts and using 
wheat Rubisco as a reference normalized to 100, following the pro- 
cedures described in Galmes et al. (2005). The Michaelis-Menten 
constants for C0 2 (K c ) and 0 2 (K 0 ) and the maximum rate for the 
carboxylase reaction (V^x) were measured in rapidly isolated leaf 
protein extracts (Sharwood et al, 2008). Briefly, ~0.5g fresh weight 
of leaves were ground in a mortar with 2 ml of ice-cold extraction 
buffer containing lOOmM Bicine (pH 8.2), 6% (w/v) polyethylene 
glycol (PEG) 4000, 2mM MgCl 2 , 0.1 mM EDTA, ImM benzami- 
dine, 1 mM e-aminocaproic acid, 50 mM 2-mercaptoethanol, lOmM 
dithiothreitol (DTT), 2 uM pepstain A, 10 uM E64, 10 uM chy- 
mostatin, 2mM phenylmethylsulphonyl fluoride (PMSF), and 2.5% 
(w/v) polyvinylpolypyrrolidone (PVPP). The liquidized sample 
was clarified by centrifugation at 12 842 g for 4min. A 1 ml aliquot 
of clarified extract was eluted using a Sephadex PD-10 column 



(GE Healthcare, UK) pre-equilibrated with desalt buffer contain- 
ing lOOmM Bicine (pH 8.2), 20mM MgCl 2 , lOmM DTT, ImM 
KH 2 P;, 0.5 mM EDTA, ImM benzamidine, ImM E-aminocaproic 
acid, and lOmM NaHC0 3 . The protein peak (in 1 ml) was supple- 
mented with protease inhibitors (4 uM pepstain A, 20 uM E64, and 
20 uM chymostatin). Of this extract, 250 ul was supplemented with 
2.542 uCi NaH 14 C0 3 and activated for 15min before carboxylase 
measurements. The remainder was used to assay Rubisco catalytic 
sites. Measurements of K c and K cmax for the carboxylase activity in 
nitrogen or air were determined in two sets of eight vials from the 
amount of 14 C incorporated into PGA, as described elsewhere (Bird 
et al, 1982). Each set of vials used eight different concentrations of 
bicarbonate chosen to provide C0 2 (aq) between 0.7 uM and 75 |iM, 
each with a specific radioactivity of 3.7 x 10 10 Bq mol -1 and con- 
taining 375 nmol RuBP. K a was calculated from the relationship K c 
(air)=K c (N 2 )x{l+[0 2 ]AK 0 }. The concentration of Rubisco catalytic 
sites in the extract was measured from the stoichiometric binding 
of the inhibitor [ 14 C]CABP to C0 2 -Mg 2+ -activated Rubisco active 
sites (Butz and Sharkey, 1989). Thereafter, the carboxylase catalytic 
turnover rate K CM C was obtained as K csl - F cmax /[catalytic sites]. The 
Rubisco kinetic constants are summarized in Table 1. 

Statistical analyses 

Data are expressed as the means ±standard error. Student's r-tests 
were used to compare the photosynthetic parameters calculated 
with the different sets of Rubisco kinetic constants and to examine 
whether the intercepts of the regression were significantly different 
from zero. Linear regression and statistical analyses were carried out 
using Microsoft Excel. 

Results 

The effect of the different calibration methods on a|3 

The four calibration methods tested here comprised scenarios 
ranging from a low (Yin PPFD <400 umol nr 2 s~') to a high 
reductant supply (O PSII /<i> C02 PPFD >400 umol nr 2 s~'), all 
performed under non-photorespiratory conditions, although 
potentially affected by the occurrence of an alternative elec- 
tron flow. Additionally, three species covering a large range in 
SLA (33 ±1.4, 14 ±2.7, and 9.0 ± 1.3 m 2 kg" 1 for N. tabacum, 
C. arabica, and L. gibertii, respectively) were selected to assess 
the extent to which structural changes might play a role in 
determining the a(3 product. 

High regression coefficients were obtained for all of 
the relationships based on the four calibration methods 
(r 2 =0.9 1-0.98) (Fig. 2). In contrast to expectations, no major 
differences in the estimation of ap were found between the 
methods Yin (PPFD <400 umol nr 2 s 1 ) and 0 PSII /O C o 2 
(PPFD >400 umol nr 2 s 1 ) (Fig. 2B, C) and between O psii / 
0 C02 (A/Q) and Yin (PPFD <400/C; >500 umol mol 1 air) 
(Fig. 2A, D). Thus, the inclusion of data at high PPFD or low 
Q, conditions that favour the prevalence of alternative elec- 
tron fluxes, had slight effects on the estimated a(3. Instead, 
the major differences among the calibration methods were 
dependent on the type of data used to perform the calibra- 
tion [i.e. ^4/PPFD or A/Q curves (Table 2)]. Regardless of the 
studied species, the calibrations using only /1/PPFD data gave 
lower a|3 values, ranging from 0.37 to 0.50, whereas those 
using A/Q data produced higher a|3 values, ranging from 0.46 
to 0.62 (Table 2). Notably, the species displaying the most 
contrasting SLA values (N. tabacum and L. gibertii) showed 
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Fig. 1 . Graphical representation of the method utilized to retrieve the respiration value from the combined A/Q and A/C c curves 
(R A Ci/Ac c )■ The example is based on values measured in an N. tabacum leaf. (A) The initial part of an A/Q (filled circles) and the 
respective A/C c (open circles) curve considering R L =R da M2 0 -1 Rmol C0 2 rrf 2 s" 1 ). (B) Q-C c from the curves presented in (A) plotted 
against A. They-axis positive intercept (2.4) when x=0 (see the dotted lines) means negative g m tor A lower than the intercept because 
Q >C C . The slope of the linear regression (0.276) is an estimation of averaged g m over the range of Q used in the linear fit. (C) By 
adjusting to zero the intercept in the relationship A versus Q-C c , a new respiration value, Raq/ac c > is obtained (0.4 umol C0 2 nr 2 s~ 1 ), 
and the same C0 2 compensation point is now shared by both curves (C-C c =0 at ^=0). The influence of the new respiration on g m can 
be observed as the modified slope (0.259). (D) The modified initial part of the/VC, and^\/C c curves calculated with R ac ,/ac c ■ 



similar values for a(3, suggesting that leaf structure does not 
play a major role in determining a(3. 

Mesophyll conductance as affected by the J 
calibration methods and respiration estimations 

Given that no major differences were found between the 
0p SII /O CO2 relationship and the Yin method, g m was esti- 
mated using the methods covering the extremes of a(3 values, 
namely O PSII /<I> C02 (A/Cd and Op SII /O C02 (PPFD >400 umol 
nr 2 s~'). Because these methods consider only ^4/PPFD 
(O PSII /O C02 , PPFD >400 umol nr 2 s" 1 ) or A/Q data [O psii / 



®co 2 the g m estimated using either method is hereaf- 

ter referred to as the g m obtained with the ^/PPFD or A/Q 
J calibration. Once the two / calibrations to be used were 
defined, an R L value next had to be chosen to calculate C c . 
Because there was no consistency in the R L estimated via the 
Yin approach using ^4/PPFD (PPFD <400 umol nr 2 s" 1 ) or 
A/Q (PPFD <400 umol nr 2 s" 1 or Q >500 umol moF 1 air) 
(Table 2), i? dark / 2 was preferred because it is the unique respi- 
ration value actually measured in planta. To check the con- 
sistency of the respiration value used, a new approach was 
also tested to find an alternative proxy for R L (Raqiac c ) as 
the value forcing the C0 2 compensation point (T) to be equal 
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Table 1 . Rubisco kinetic constants measured for the species studied: 
specificity factor (S c/a ), C0 2 compensation point in the absence 
of respiration (T*), Michaelis-Menten kinetics for C0 2 (KJ and 0 2 
(KJ, and catalytic turnover rate for the carboxylase reaction (K cat ) 



Species 


S c /o 


r* (nbar) 


K c ( M M) 


K 0 ( M M) 


K cat (s- 1 ) 


N. tabacum 
C. arabica 
L. gibertii 


98.1 + 2.6 a 
98.4 + 4.3 
11 0.5 ±1.6" 


39.7 + 1.1 
39.6+1.7 
35.2 + 0.5 


12.4 + 0.7 
10.3+1.3 
8.9 + 0.5 


274 + 42 
479 + 113 
593 + 75 


3.2 + 0.3 
3.2 + 0.1 
2.7 + 0.8 



Values are the means ±standard error of 3-4 replicates per species. 
a Taken from Galmes et al. (2006) 
b Taken from Galmes et al. (2005). 



when calculated from A/Q and A/Q curves. The estimated 
R A qiac c values were always lower than their corresponding 
^dark/2 counterparts (Table 3). 

Due to the high amount of negative and extremely high g m 
values at low or high Q, respectively (Table 4), a filter was next 
applied to keep the valid g m estimates, here defined as those val- 
ues in the range of 0 <g m <1 mol C0 2 nr 2 s" 1 . As expected, the 
use of Raqiac c significantly reduced the amount of negative 
values of g m at low Q (Raqiac c i n Table 4), as the rationale 
for this method works at this Q range (Fig. 1). Conversely, the 
use of A/Q / calibration improved the g m estimation at high Q 
(Table 4) due to a lowering of Q values (Fig. 3A). Thus, the 




PPFD4 psll /4 PPFD* psll /4 



Fig. 2. Calibration relationships O PS h versus <£ C o 2 (A and B) and A versus PPFD fl>p SN /4 (C and D) measured under non- 
photorespiratory conditions (<1 % 0 2 ) by varying PPFD intensities (/4/PPFD curves) or substomatal C0 2 concentrations (/A/C, curves). For 
more details on the four calibration methods, see the Materials and methods. In (A), the entire A/C t curve is utilized; in (B), only values at 
PPFD >400 ixmol rrr 2 s~ 1 are considered; in (C), only values at PPFD <400 |imol nr 2 s~ 1 are used; in (D), only values at PPFD <400 nmol 
rrr 2 s _1 plus C, >500 nmol mol" 1 are considered. The slope of the lines in all graphs refers to the product a|3, whereas the y-intercept 
should be interpreted as the presence of alternative electron sinks in A and B and as a measure of R L in C and D. The values of the 
slopes and intercepts are summarized in Table 2. Only <1^ 02 values <0.05 were kept in the plots, as recommended by Seaton and 
Walker (1990). The scales have different amplitudes according to the calibration. In the plots, data from 4-6/VPPFD orA/C; curves are 
used for the linear regression. 
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Table 2. Slopes (a/3) and light respiration following the method of Yin et al. (2011) (R L ym, /umol C0 2 rrf 2 s' 1 ) or intercept values (0 PS i/ 
%q 2 ) obtained under non-photorespiratory conditions according to different approaches, consisting of higher (<Z> PS // <% 02 A/C, or 
PPFD >400 fjmol rrr 2 s~ 1 ) or lower (Yin PPFD <400 or PPFD <400 fjmol m~ 2 s~ 7 and C, >500 /jmol mol~ 1 air) susceptibilities to alternative 
electron sinks 

For comparison, the original slopes (k) obtained when using the Opsi/^bc^ relationship were already converted to a(3 (a(3=4//c); the intercept 
refers to the parameter b in the equation: <D PS u = k 4> C02 +b. 







Slope (ap) 


Intercept (a>p SM / <t> C02 ) 




N. tabacum 


^PSIK ^C02 V" 


0.52 + 0.013 


0.005 + 0.002 








0.36 + 0.011 


-0.02 + 0.008* 






Yin (PPFD <400) 


0.39 + 0.017 




0.87 + 0.47 




Yin (PPFD <400 C, >500) 


0.46 + 0.026 




0.87 + 0.87 


C. arabica 


1>psii/<1>co 2 WQ 


0.64 + 0.014 


0.015 + 0.002* 






<W *co 2 (PPFD>400) 


0.46 + 0.018 


-0.011+0.007 






Yin (PPFD <400) 


0.50 + 0.025 




0.84 + 0.41 




Yin (PPFD <400 C, >500) 


0.63 + 0.015 




2.55 + 0.45* 


L. gibertii 


fpSl/ < I > C022 (/A/C|) 


0.56 + 0.019 


0.014 + 0.004* 






1>psii/<1>co 2 (PPFD>400) 


0.37 + 0.011 


-0.03 + 0.009* 






Yin (PPFD<400) 


0.41 +0.013 




2.04 + 0.25* 




Yin (PPFD <400 C, >500) 


0.50 + 0.018 




3.14 + 1.02* 



Values are the means +standard error of four replicates per species. 

An asterisk denotes respiration or intercepts significantly different from zero (P < 0.05). 



Table 3. Dark respiration measured at pre-dawn (R dark ) and light 
respiration estimated from combined A/C, and A/C c curves 

f^ACj/ACc ) 

All values are in ^mol C0 2 nr 2 s~ 1 . 



Species 


Rdark 




N. tabacum 


2.2 + 0.1 


0.3 + 0.1 


C. arabica 


0.9 + 0.1 


0.2 + 0.1 


L. gibertii 


4.0 + 0.3 


1.7 + 0.4 



Values are the means ±standard error of four replicates per species. 
a For R AC/AC , the values are averages from the A/C c curves 

calibrated using the,4/PPFD or A/C-, curves. 

highest percentage of data excluded (i.e. unrealistic g m values) 
was found when combining the A/PPFD J calibration with 
-^dark/2 (88, 56, and 42% for N. tabacum, C. arabica, and L. giber- 
tii, respectively), whereas the lowest percentage of data excluded 
was obtained with the Al QJ calibration and Raqiac c (17% for 
N. tabacum and 0% for the other species) (Table 4). 

The different calibrations also affected the magnitude 
of g m in response to Q. There was no general trend at Q 
<100 umol mol" 1 air, whereas at Q >100 ■xmol mo\ [ air, 
lower g m values (up to 50%) were observed using the A/Q J 
calibration (Table 4). In addition, when plotting g m versus Q, 
a high amplitude in g m , depending on the calibration used, 
was revealed (Fig. 3B). Such amplitude led to a reasonable 
disagreement between the averaged g m at Q of 100-350 umol 
moh 1 air and the single point g m at C0 2 ambient (400 umol 
moh 1 ) (Supplementary Tables S2, S3 at JXB online), which 
suggests that care must be exercised when reporting single 
point g m data. In Supplementary Table S4, in addition to the 



filter keeping the values in the range of 0 <g m < 1 mol C0 2 
rrr 2 s" 1 , the Harley et al. (1992) criteria were also applied, 
considering the g m data in the range of 10 <dCJdA <50 as 
reliable to see whether they corresponded to the data pre- 
sented in Table 4. Overall, the g m values presented in Table 4 
and those obtained after applying the Harley criteria varied 
accordingly, but some exceptions were noticed such as the 
g m values estimated at C, >350 umol mor 1 air for L. giber- 
tii which were substantially higher than the g m data shown 
in Supplementary Table S4. Most importantly and irrespec- 
tive of this difference, the amount of data excluded using the 
Harley criteria was always higher (up to 100%) even when the 
same estimates presented in Table 4 had 0% of data excluded 
together with small standard errors. Additionally, the Harley 
criteria worked differently depending on the species, / cali- 
bration, and Q range, given that some species were more 
affected than others. 

To obtain an independent estimate of g m that was not 
affected by the need for calibrating /, g m was also estimated 
using an alternative approach, namely the A-Q curve analy- 
sis method suggested by Ethier and Livingston (2004). The g m 
values obtained following this approach are given in Table 5. 

Maximum Rubisco carboxylation and electron 
transport rate as affected by a|3 

In contrast to the g m estimates that could be greatly affected 
by the use of either or R ACi i A c c (Table 4, Fig. 3B), 

the F cmax and / max estimates were minimally affected by the 
use of different respiration values (~6% at most; Table 5; 
Supplementary Table S5 at JXB online). 

Irrespective of species, F cmax and / max on a C c basis were 
always higher for the A/Q / calibration than for its ^4/PPFD / 
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Table 4. Mesophyll conductance (g m mol C0 2 m 2 s 1 ) for several intervals of C, and percentage of data excluded (DE) after applying a 
restriction (g m restricted to the range ofO <g m <1 mol C0 2 rrf 2 s~ 1 ) 

TheA/Q or/VPPFD J calibration refers to the Opsi/^boo (4/Q) or/4/PPFD methods, respectively. /-?aq/ac c ' s the ''9 nt respiration estimate from the 
combined A/C, and A/C c curves proposed in this study, and R dar k/2 is the dark respiration divided per two to account for the observed reduction in fl dark 
under light (Niinemets ef al. , 2005). Note the high DE at low Q for R dar k/2 and at high Q for the A/PPFD J calibration. 





A/C, J calibration 








yVPPFD J cal 


ibration 






"ACi/ACc 




"dark/2 




^ACi/ACo 




"dark/2 




9m 


DE \/o) 


dm 


r\c to/ \ 
DE (%) 


0m 


r\c 10/ \ 
DE (%) 


dm 


r\c 10/ \ 
DE (%) 


N. tabacum 


















Ci<100 


0.290 + 0.065 


17 


0.243 + 0.033 


58 


0.345 + 0.132 


17 


0.203 + 0.031 


67 


Q 100-350 


0.213 + 0.021 


0 


0.287 + 0.037 


6 


0.468 + 0.064 


50 


0.504 + 0.262 


88 


Q >350 


0.038 + 0.005 


0 


0.045 + 0.007 


0 


0.262 + 0.024 


65 


0.132 + 0.130 


70 


C. arabica 


















Ci<100 


0.131 +0.018 


0 


0.202 + 0.042 


10 


0.164 + 0.028 


0 


0.202 + 0.035 


20 


Q 100-350 


0.117 + 0.006 


0 


0.124 + 0.007 


0 


0.186 + 0.029 


0 


0.191 +0.029 


6 


Q >350 


0.052 + 0.004 


0 


0.053 + 0.004 


0 


0.219 + 0.038 


41 


0.194 + 0.039 


56 


L gibertii 


















Ci<100 


0.235 + 0.027 


0 


0.270 + 0.050 


7 


0.310 + 0.041 


7 


0.275 + 0.079 


27 


Q 100-350 


0.214 + 0.014 


0 


0.228 + 0.020 


0 


0.371 +0.044 


7 


0.406 + 0.083 


21 


Q >350 


0.071 +0.011 


0 


0.073 + 0.012 


0 


0.226 + 0.063 


32 


0.162 + 0.044 


42 



Values are the means ±standard error of 4-6 A/Q curves per species. The SE was calculated according to the points that remained in each Q 
interval after applying the restriction. 



counterpart (Table 5). These higher F cmax and / max values are 
in agreement with the lower g m values found, in general, with 
the Al Q / calibration (Table 4). The F cmax obtained using the 
Ethier and Livingston method was closer to the F cmax on a 
Q basis calculated using the Al Q J calibration for L. gibertii 
and C. arabica, whereas for N. tabacum an intermediate value 
between those obtained with the ^4/PPFD or AIQ J calibra- 
tion was found (Table 5). Apart from the / max calculated with 
the A/Q / calibration, which had the higher values for all spe- 
cies and reflected the higher ap found (Table 3), there were no 
major differences among the other calculated values of / max . 
The use of the standard Rubisco kinetics (K c , K 0 , and T*) 
originally obtained for N. tabacum by Bernacchi et al. (2002) 
would lead to an overestimation of Fc max by ~30% and 20% 
for L. gibertii and C. arabica, respectively. The other photo- 
synthetic parameters for these two species were unaffected 
by the use of different Rubisco kinetic constants (Table 5). 
No significant differences were observed for N. tabacum when 
using the different set of Rubisco kinetics, with the exception 
of g m (Ethier and Livingston method), which was 28% lower 
when using the Rubisco kinetics of Bernacchi et al. (2002). 



Discussion 

To the best of the authors' knowledge, this study is the first to 
address the idea that the use of A/Q or A/PFFD curves under 
low 0 2 to calibrate / can significantly affect the estimations of 
gm, ^cmax* and / max , as demonstrated in species with contrasting 
leaf structural properties and photosynthetic capacities. A back- 
ground is also provided to understand potential factors that 
could be translated into unrealistic g m values at low and high Q, 
and recommendations to improve g m estimations accordingly. 



Sensitivity of the variable J method to a|3 

The high sensitivity of the variable / method has been known 
since its introduction by Harley et al. (1992). However, 
among the several sources of error, more attention has been 
given to an accurate determination of T* (Warren et al., 
2006; Pons et al. , 2009) than to / per se, given the variety of 
ways in which / can be calibrated. Despite the / calibration 
based on ^4/PPFD curves being the most common method 
in the literature, here it is shown that this calibration pro- 
duced a higher number of unrealistic g m estimates than its 
counterpart based on A/Q curves (Table 4). Additionally, it 
was demonstrated that g m misestimations could be directly 
associated with the / calibration method at high Q on the 
one hand, and the proper choice of T* and R L at low Q on 
the other hand. Notably, the need to develop useful criteria 
to improve g m estimations and to understand the high g m 
variability is crucial considering that the only existing indi- 
cator, the Harley et al. (1992) criteria which consider that 
reliable g m estimations are those situated in the range of 
dC c /d^4 (10-50), performed poorly in retrieving acceptable 
data in addition to being dependent on the species, applied 
Q range, and / calibration (Supplementary Table S4 at JXB 
online). 

It can be deduced from the relationship between Q and 
J/(A+R L ) (Fig. 3A) that Q approaches infinity as J/(A+R L ) 
approaches four (number of electrons per C0 2 molecule 
fixed). In contrast, Q is negative when the ratio is less than 
four [the lower Q threshold for positive J/(A+R L ) is defined 
by the parameter T*]. Therefore, as J/(A+R L ) approaches 
four as photorespiration and other alternative electron sinks 
tend to decrease (at high Q, for example), extremely high C c 
values are to be expected, leading to a greater probability of 
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Fig. 3. (A) Relationship between C c and J/iA+Rf), as affected by the J calibration using the A/PPFD (open circles) or A/C t curves (filled 
circles). As J/IA+R^ approaches four (likely to occur at high Q), C c tends to infinity, thus increasing the probability of being higher than 
Q and resulting in negative g m values. The arrow indicates the C c upper limit (-600 nmol mol" 1 ) observed for the/A/C, J calibration, 
being much lower than its A/PPFD counterpart which can reach C c values >1 600 jimol mol" 1 . (B) g m response to C, as affected by the 
/A/Cj (filled symbols) or A/PPFD (open symbols) J calibration method and the light respiration estimations F? dark / 2 (triangles) or R ACVAC 



(circles)]. It is remarkable how there are 'spikes' in g m at low Q when using R da M2 and how these are alleviated when using ft 

curves are from L. gibertii, and the points are averages from four plants for g m between 0 and 1 mol C0 2 rrr 2 s~ 
calibration refers to the a|3 obtained from the curves shown in Fig. 2A and B, respectively. 
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Table 5. Maximum carboxylation rate ofRubisco (V cmax , pmol C0 2 m~ 2 s~'J, maximum electron transport rate from gas exchange 
(J max , fjmol e~ m" 2 s~ 1 ), mesophyll conductance (g m Ethier, mol C0 2 m~ 2 s" 1 ) according to Ethier and Livingston (2004), and electron 
transport rate estimated from chlorophyll fluorescence (Jflu, /jmol e~ m~ 2 s~ 1 ) 

All of the photosynthetic parameters were calculated on a C, basis from the A/Q curves or on a C c basis from the A/C c curves using the 
measured Rubisco kinetic constants (T*, K c , and K 0 ) in this study (as reported in Table 1) or the standard Rubisco kinetics for tobacco from 
Bernacchi ef al. (2002). A/Q J cal and /4/PPFD J calibration denote the photosynthetic parameters calculated for A/C c curves using the af> 
retrieved from <D PSII / o COz (4/Cj) or (4/PPFD) J calibrations. Note that J flu is independent of Rubisco kinetics. 



Measured kinetic constants Bernacchi constants 





N. tabacum 


C. arabica 


L. gibertii 


N. tabacum 


C. arabica 


L. gibertii 


^cmax Q basis 


71.7 + 4.6 


41.7 + 2.4* 


70.5 + 1.6* 


59.8 + 3.8 


49.8 + 2.7 


104.0 + 2.4 


Kxnax C c basis >VCi J Cal 


106.3 + 8.9 


73.9 + 3.6* 


133.4+11.9* 


87.0 + 7.1 


89.8 + 4.4 


210.1 +20.8 


^ ma xC cbasis /\/PPFDJcal 


80.2 + 5.7 


59.8 + 3.8* 


101.7 + 8.2* 


66.7 + 4.6 


70.3 + 3.1 


155.8+14.0 


V cmax Ethier 


91.8 + 14.7 


80.2 + 4.0* 


124.8 + 5.0* 


81.7 + 9.2 


102.1 +5.4 


164.6+11.7 


^max basis 


87.2 + 7.4 


110.7 + 3.4 


175.3 + 7.3 


86.6 + 7.4 


110.7 + 3.4 


176.6 + 7.4 


^max C c basis A/C\ Cal 


103.2 + 10.5 


136.6 + 7.1 


228.6 + 8.7 


101.4 + 10.3 


136.6 + 7.1 


233.4 + 8.9 


^max C c basis /A/PPFD Cal 


84.9 + 7.3 


113.7 + 5.7 


177.0 + 5.8 


84.4 + 7.3 


113.1 +5.3 


178.2 + 5.8 


Jmax Ethier 


99.5 + 3.3 


115.0 + 3.2 


179.8 + 8.9 


99.5 + 3.3 


114.4 + 3.2 


179.8 + 8.9 


J, u A/C t J cal 


117.5 + 11.4 


139.2 + 6.7 


236.0 + 9.7 








J flu /VPPFD J cal 


94.7 + 7.9 


113.1 +4.8 


181.7 + 6.4 








9m Ethier 


0.412 + 0.045 


0.108 + 0.006 


0.226 + 0.036 


0.296 + 0.015 


0.101 +0.005 


0.257 + 0.048 



Values are the means ±standard error of 4-6 A/C, orA/C c curves per species. 
Raciiac c was used when needed. 

An asterisk denotes differences in the respective photosynthetic parameter calculated using the different set of Rubisco kinetics. 
The graphical representation of the A/C, and A/C c curves is shown in Supplementray Figure S1 at JXB online. 



estimating negative values of g m . Thus, for the same (A+R L ), 
the higher the a|3, the higher the /, ultimately resulting in a 
greater probability of obtaining positive g m values due to a 
lowering of Q, whereas the opposite holds true for a lower 

At low Q (<100 umol mol" 1 air), misestimations of g m 
are unlikely to be affected by J/(A+R L ), as this ratio tends 
to achieve higher values with decreasing Q. In any case, reli- 
able estimations of g m at low Q are extremely challenging 
given that the estimations are highly dependent on the proper 
choice of values of T* or R h , which can ultimately affect C c . 
In the present study, the authors are quite confident about 
their F* values because they were estimated from Rubisco 
kinetics in purified preparations. Given this, it was possible to 
retrieve a respiration estimate (R A Cj/Ac c ) as the value forcing 
the intercept of the linear relationship of A versus Q-Q to 
zero, which is equivalent to forcing the F to be the same. It 
is believed that R ACi /Ac c obtained in this way is a valid esti- 
mate, as its magnitude was lower than that of R dMk (Table 3), 
in agreement with the general consensus that respiration is 
inhibited in the light (Tcherkez et al., 2005). In addition, the 
use of RaqiaCc significantly improved the number of valid 
g m estimates (Table 4), highlighting the importance of respi- 
ration in estimating g m at low Q, with the advantage of bet- 
ter matching a biochemical prediction of the FvCB model, 
namely A/Q and A/C c sharing the same T. 

An a priori weakness regarding Raq/aCc resides in whether 
g m is kept constant at a low Q range. However, the authors 
argue in favour of a constant g m where a high linearity can be 
observed when plotting A versus Q-Q, in the region strictly 
limited by Rubisco (Fig. 1C). Furthermore, it can also be 



noted that g m estimated using the variable / method is closer 
to the constant g m , as the individual points are close to the 
regression line in the plot of A versus C—C c . The apparent 
variability of g m that might be identified using the variable / 
method could simply be a result of the high sensitivity of A, 
Q, or C c to random errors at this low Q range. This conclu- 
sion is supported when comparing the data of Flexas et al. 
(2007), who showed a high variability of g m at low Q using 
the variable / method, with those recorded by Tazoe et al. 
(201 1 ), who found almost constant values of g m using the iso- 
topic method. 

Why do A/C, and A/PPFD curves result in different a|3? 

Given that a can be measured and that large changes in this 
parameter during the execution of A/PPFD or A/Q curves 
are unlikely to occur, the uncertainty in / is particularly 
related to |3. Eichelmann and Laisk (2000) found (3 vary- 
ing from 0.38 to 0.51 in N. tabacum under varying light and 
temperature conditions. Similar results were reported by 
Loreto et al. (2009) while studying the effect of blue light 
on g m . Hassiotou et al. (2009) also concluded that / calibra- 
tion is light dependent. Collectively, all of this information 
highlights the importance of irradiance in determining the p 
value and helps explain the apparently better suitability of 
the .4/C r based / calibration: it is performed under the same 
irradiance used in the normal A/Q curves. In fact, another 
advantage of keeping a fixed irradiance during the calibration 
is to avoid changes in the profiles of light absorption through 
the leaf that may occur under changing irradiance (Evans, 
2009; Oguchi et al. , 201 1); this can be especially important for 
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the Li-Cor LED-based fluorescence/light source which con- 
sists of blue and red light with narrow bandwidths. However, 
to what extent changes in spectral light could be responsible 
for differences in |3 is unclear (Evans, 2009). Further explana- 
tions for this statement might be linked to the engagement 
of alternative electron sinks and the fluorescence signals that 
might not be representative of the whole leaf, in contrast to 
the gas-exchange signals (Hassiotou et al, 2009). Even if the 
uncertainty about the fluorescence signal could be resolved 
using calibration curves, as in this study, conflicting informa- 
tion on the possible effects of alternative electron sinks on 
g m measurements has been reported. Whereas several authors 
(e.g. Loreto et al, 1994; Ruuska et al, 2000; Flexas and 
Medrano, 2002) showed circumstantial evidence suggesting 
that alternative electron sinks play no major role in g m esti- 
mations, other investigators reported that up to 24% of the 
total electron flux can be associated with alternative sinks (see 
Gilbert et al, 2012, and references therein). In any case, the 
present results argue against the relevance of alternative elec- 
tron sinks as potential bias for proper g m estimations, given 
that the calibration methods most likely to be affected by 
these sinks (calibrations performed at high PPFD and low Q) 
did not differ from others which used a range of data under 
strictly limited electron transport (low PPFD and high Q). 

Mesophyll conductance, maximum velocity of Rubisco 
carboxylation, and electron transport rate 

The choice of the method to calibrate / can significantly 
affect the estimation of photosynthetic parameters and even 
alter their interpretation in the context of the diffusive versus 
biochemical limitations to photosynthesis. In N. tabacum, the 
use of ^4/PPFD calibration would lead to the conclusion that 
the assumption of infinite g m is plausible because no differ- 
ence in F cmax on a Q or C c basis was observed (Table 5). In 
sharp contrast, the use of Al Q calibration points to a finite 
g m and a higher F cmax value (33%) in relation to the F cmax on 
a Q basis. For C. arabica and L. gibertii, higher F cmax values 
were observed using both A/Q and ^4/PPFD / calibrations 
compared with F cmax on a Q basis, but the degree of under- 
estimation varied considerably (Table 5). It is important to 
bear in mind that such alterations also imply changes in bio- 
chemical aspects of the leaf because F cmax is related to the 
amount of activated Rubisco, whereas / max is associated with 
components of the electron transport chain, and both param- 
eters can affect the leaf nitrogen economy (Niinemets and 
Tenhunen, 1997). Thus, it is highly unlikely that both calibra- 
tions can hold true, as they can considerably change the leaf 
biochemical signature and hence ultimately pose uncertain- 
ties on how to decide which calibration to use. To address this 
issue, a first cut-off point is to check if the minimum require- 
ment of four electrons per carboxylation is attained, which 
can be achieved by dividing / by A. If this ratio is less than 
four, the calibration is, on a theoretical basis, inadequate. 
This criterion allowed the exclusion of the standard a(3 for 
C. arabica (Supplementary Table S2 at JXB online) and the 
A/PPFD calibration for N. tabacum. In addition, g m and the 
FvCB photosynthetic parameters were estimated using the 



/-independent Ethier and Livingston method. Interestingly, 
the values of F cmax and g m obtained with this method rea- 
sonably matched those obtained with the A/Q J calibration 
for L. gibertii and C. arabica. In N. tabacum, the F cmax calcu- 
lated from the Ethier approach presented intermediate values 
between the F cmax from the ^/PPFD or A/Q J calibration, 
and, given the high standard errors in its g m estimate, it is 
believed that, even in this species, the A/Q / calibration might 
be a better option to recommend, as it allows the retrieval of 
a greater number of realistic g m estimations regardless of the 
respiration value (Tables 4, 5). Therefore, it is proposed that 
the A/Q J calibration seems to be more reliable than its Al 
PPFD counterpart in terms of producing acceptable data, in 
accordance with the results of Gilbert et al. (2012). 

Regarding the use of Rubisco kinetic constants, a consider- 
able overestimation of F cmax (~30%) would occur in /. gibertii 
if standard values of those constants (Bernacchi et al. , 2002) 
were used (Table 5). Importantly, even if good fits could be 
obtained irrespective of the kinetic constants, the relation- 
ship between V cm!lx and leaf nitrogen might be compromised. 
Thus, if the main goal is to use the FvCB model to charac- 
terize photosynthetic capacities and relate them to nitrogen 
partitioning (e.g. Xu et al, 2012), species-specific kinetic con- 
stants for Rubisco should be implicitly used. 

Further recommendations to improve g m estimation 

Given the uncertainties in g m estimations, criteria were pro- 
vided to assess the consistency of T* and R L through the 
determination of T (by definition, T must be higher than 
T*) or by plotting A versus Q-C c and analysing the inter- 
cept of the linear relationship. In addition, provided that the 
estimates of T and T* are reliable, it is possible to retrieve a 
respiration estimate from A/Q and A/C c curves. The respira- 
tion value was the focus here because of the availability of 
the current methods to estimate R L using gas exchange and/ 
or chlorophyll fluorescence (see Yin et al. , 201 1), all of which 
were performed under low irradiance conditions that do not 
match those used in A/Ci curves, which ultimately makes the 
choice of a proper respiration value a very complicated task. 

If an inconsistency between T for the A/Q curve and T* 
(T lower than T*) is found, or if the value of respiration 
retrieved by forcing the intercept of A versus C—Q to zero is 
negative, four potential sources of error need to be checked: 
(i) correction for C0 2 and water vapour leakage through the 
chamber gaskets is crucial to the determination of T and 
very important for species with low photosynthetic poten- 
tial (see Rodeghiero et al, 2007); (ii) the influence of lateral 
diffusion which can become significant (especially for homo- 
baric leaves) with large gradients of C0 2 inside and outside 
the chamber in addition to the lower fluxes of C0 2 near the 
compensation point (Morison and Lawson, 2007); (iii) T* 
determination; and (iv) Rubisco deactivation. The first two 
sources of error are particularly dependent on the leaf cham- 
ber size and, given the magnitude of the corrections in A at 
low Q (Supplementary Table SI at JXB online), the need to 
use larger leaf chambers to obtain more reliable g m , T, and 
RaqiaCc estimates is emphasized. For species displaying low 
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respiration rates, it is possible to obtain an estimate of T* 
by fixing R L , inputting a test T*, and changing this param- 
eter to obtain the intercept of A versus C~C C closest to zero. 
All of these hints provide a good framework to extract as 
much information as possible from the A/Q and A/Q curves. 
Nevertheless, because all of these adjustments rely on fitting 
processes, proper judgement about the biological meaning of 
the obtained estimates is obviously required. 

Regarding the / calibration, independent of it being based 
on A/Q, ^4/PPFD, or both, one should first check whether all 
of the positive JIA values are higher than four. If not, there is 
an inconsistency with the calibration because this implies that 
fewer than four electrons would be used per carboxylation 
event. If yes, the next step is to verify the suitability of the res- 
piration value by checking if all of the positive J/(A+R L ) val- 
ues are still higher than four. If not, there is an inconsistency 
for the same reason stated previously. Another checkpoint is 
to verify that the maximum observed Q is at least equal to 
the maximum Q that would implicate infinite g m . If C c values 
higher than the maximum Q are present for positive A, there 
is again an inconsistency with the calibration, as Fick's first 
law of diffusion predicts a drawdown of C0 2 between Q and 
Q. In fact, it is theoretically possible to define an upper limit 
for ct|3 which would be the value making the Al Q and Al Q 
curves identical, implicating infinite g m . 

Conclusion 

This work brings new information to the growing amount of 
published papers trying to improve g m estimation. Criteria 
are provided to examine the amount of unrealistic g m data 
and identify inconsistencies with the biochemical model 
utilized, leading to a higher amount of acceptable data. In 
addition, the need for additional measurements to validate 
the / calibration to improve the estimation of photosynthetic 
parameters is emphasized. Moreover, given the observed high 
variability in g m , it is important that comparisons among 
studies reporting data consider the variability due to the / 
calibration method, and, if necessary, normalize data accord- 
ing to the reported ap values. 



Supplementary data 

Supplementary data are available at JXB online. 

Figure SI. Graphical representation of the A/Q and their 
respective A/C c curves calibrated with the A/Q or A/PPFD J 
calibration. 

Table SI. Leakage correction for C0 2 and water vapour. 

Table S2. Averaged mesophyll conductance (g m ) for the 
interval of Q ranging from 100 to 350 umol mok 1 air and 
the single point g m at ambient C0 2 concentration (400 umol 
mok 1 air) using -R dark / 2 in the g m estimation. 

Table S3. The same as in Table SI, using R A qiac c rather 
than -R dai . k/2 to estimate g m . 

Table S4. Mesophyll conductance (g m , mol C0 2 nr 2 s~') for 
several intervals of Q and percentage of data excluded (DE) 
after applying two restrictions (g m restricted to the range of 



0 < g m <lmol C0 2 nr 2 s 1 and dQ/dA of 10-50 [the Harley 
et al. (1992) criteria]. 

Table S5. Maximum carboxylation rate of Rubisco (K cmax ) 
and maximum electron transport rate from gas exchange 
(•Anax)> as affected by the different / calibrations and calcu- 
lated with i? dark/2 . 

Spreadsheet SI. An Excel spreadsheet is provided that 
allows the user to retrieve the leaf respiration from combined 
A/Q and AIC C curves. 
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